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Abstract 

A  thin  cylindrical  ceramic  sample  is  placed  in  a  single  mode  microwave 
applicator  in  such  a  way  that  the  electric  field  strength  is  allowed  to  vary 
along  its  axis.  The  sample  can  either  be  a  single  rod  or  two  rods  butted 
together.  We  present  a  simple  mathematical  model  which  describes  the  mi¬ 
crowave  heating  process.  It  is  built  on  the  assumption  that  the  Biot  number 
of  the  material  is  small,  and  that  the  electric  field  is  known  and  uniform 
throughout  the  cylinder’s  cross-section.  The  model  takes  the  form  of  a  non¬ 
linear  parabolic  equation  of  reaction-diffusion  type,  with  a  spatially  varying 
reaction  term  that  corresponds  to  the  spatial  variation  of  the  electromagnetic 
field  strength  in  the  waveguide.  The  equation  is  analyzed  and  a  solution  is 
found  which  develops  a  hot  spot  near  the  center  of  the  cylindrical  sample 
and  which  then  propagates  outwards  until  it  stabilizes.  The  propagation  and 
stabilization  phenomenon  concentrates  the  microwave  energy  in  a  localized 
region  about  the  center  where  elevated  temperatures  may  be  desirable. 
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1.  Introduction. 


The  use  of  microwaves  to  sinter  or  join  ceramics  is  rapidly  gaining  accep¬ 
tance  in  industry  where  the  efficient  production  of  high  quality  materials  is 
important.  Efficiencies  are  increased  because  microwaves  penetrate  a  ma¬ 
terial  and  rapidly  deposit  energy  there,  in  direct  contrast  to  conventional 
heating  schemes  where  heat  diffuses  into  a  material  from  its  surface.  The 
price  paid  is  the  need  of  control  systems  to  prevent  thermal  runaway  and 
other  related  instabilities. 

The  control  and  dependability  of  these  processes  requires  a  deep  under¬ 
standing  of  the  inherent  physics  which  axe  described  by  a  formidable  nonlin¬ 
ear  initial  boundary  value  problem.  This  is  comprised  of  the  time-harmonic 
version  of  Maxwell’s  equations,  the  heat  equation,  an  equation  of  state  re¬ 
lating  the  effective  electrical  conductivity  to  the  temperature,  and  a  thermal 
boundary  conditions  on  the  surface  of  the  ceramic  material  which  balances 
conduction,  convection,  and  thermal  radiation.  The  nonlinear  character 
arises  from  the  dependence  of  the  electric  field  upon  the  effective  electri¬ 
cal  conductivity,  which  is  a  function  of  the  temperature,  the  dependence  of 
the  temperature  upon  the  microwave  power  deposition,  which  is  proportional 
to  the  product  of  the  effective  electrical  conductivity  and  the  the  magnitude 
of  the  electric  field  squared,  and  the  radiative  heat  loss,  which  varies  as  the 
fourth  power  of  the  temperature.  In  addition,  the  boundary  value  charac¬ 
ter  of  the  problem  is  also  challenging  because  the  electromagnetic  fields  and 
the  ceramic  material  are  confined  in  a  cavity  or  waveguide  applicator  of  a 
complicated  geometry. 

The  systematic  analysis  of  these  equations  under  a  variety  of  physical  lim¬ 
its  has  primarily  been  restricted  to  one  dimensional  geometries  (see  reference 
4  and  the  bibliography  therein),  but  has  recently  been  extended  to  three  di¬ 
mensions  [4]  in  the  small  Biot  number  limit.  However,  in  all  these  cases 
the  effect  of  the  waveguide  applicator  or  cavity  were  neglected,  i.e.,  the  ce¬ 
ramic  samples  were  irradiated  by  plane  waves  in  free  space.  Nonetheless,  the 
small  Biot  number  theory  predicts  the  phenomenon  of  thermal  runaway  and 
suggests  methods  for  its  control.  -r 

The  problem  we  model  and  study  in  this  paper  is  concerned  with  sintering 
and  joining  of  ceramic  fibers  in  a  microwave  applicator.  As  such,  it  is  not 
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described  by  the  theories  mentioned  above,  because  of  the  applicator  and  also 
because  of  the  small  aspect  ratio  a/d  of  the  fibers.  In  this  paper,  we  take 
into  account  the  effects  of  the  applicator  by  assuming  that  the  electric  field  is 
uniform  throughout  the  cylinder’s  cross-section  and  known  along  its  length. 
That  is,  the  ceramic  cylinder  is  thin  enough  not  to  perturb  the  electric  field 
to  leading  order.  Thus,  the  heating  process  will  be  modeled  by  a  nonlinear 
heat  equation  and  boundary  condition. 

There  are  three  small  parameters  that  arise  from  a  dimensional  analysis 
of  the  simplified  model  problem.  The  first  is  the  aspect  ratio  defined  above, 
the  second  is  the  Biot  number  B\ ,  which  is  a  measure  of  convective  heat  loss 
at  the  surface,  and  the  third  is  B2  is  a  measure  of  radiative  heat  loss  there. 
As  described  above,  we  have  developed  an  asymptotic  theory  to  study  the 
microwave  heating  of  ceramic  slabs  and  other  compact  geometries  [3,4]  as 
B\  —y  0.  In  the.se  studies  B2  ~  Bi  so  that  both  physical  effects  of  radiation 
and  convection  have  been  incorporated  into  the  theory.  This  asymptotic 
theory  can  be  employed  to  analyze  the  present  problem  with  the  proviso 
that  the  parameter  =  {a/d)^/B\  is  order  one.  The  net  result  is  that 
the  temperature  remains  spatially  uniform  across  the  cereimic’s  cross-section 
and  satisfies  a  nonlinear  reaction-diffusion  equation  along  its  length.  In  this 
equation  the  reaction  term  accounts  for  adsorption  of  microwave  energy  and 
loss  of  thermal  energy,  which  arises  from  convection  and  radiation  at  the 
sample  boundaries,  and  the  diffusion  coefficient  is  c^.  The  later  will  be  taken 
as  small  (which  is  the  case  for  fibers)  to  allow  an  analysis  of  this  equation. 

Two  types  of  problems  naturally  arise  depending  upon  the  orientation 
of  the  ceramic  cylinder  in  the  applicator.  If  the  sample  is  placed  so  that 
the  electric  field  has  no  spatial  variation  Eilong  its  axis,  then  the  reaction- 
diffusion  equation  has  constant  coefficients.  Equations  of  this  type  have 
received  considerable  study  because  of  their  applicability  in  a  wide  variety  of 
physical  settings  [9].  That  the  present  equation  supports  traveling  transition 
layers  comes  as  no  mathematical  surprise.  However,  in  the  present  physical 
context  it  does  explain  the  mechanism  for  the  formation  and  propagation  of 
hot-spots  [5,6]  which  are  seen  in  experiments  [8,11]. 

On  the  other  hand,  if  the  ceramic  fiber  is  placed  so  that  the  electric  field 
varies  along  its  axis,  then  the  equation  has  a  spatially  varying  reaction  term. 
The  analysis  and  understanding  of  the  solutions  of  these  types  of  equations 
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are  less  well  understood  than  those  described  above  [l].  We  shall  analyze  the 
equation  below  using  standard  matched  asymptotic  methods  and  show  how 
its  solution  evolves  into  a  hot  siJot,  propagates  outward  from  its  inception, 
and  stabilizes  to  form  a  region  of  elevated  temperature.  This  stable  region  of 
elevated  temperature  can  be  exploited  in  fiber  sintering  and  has  already  been 
used  in  joining  processes  [2].  Moreover,  our  analysis  shows  how  the  size  of 
the  spot  depends  upon  the  temperature  dependence  of  the  material’s  thermal 
properties.  These  dependencies  are  often  ignored  because  the  experiments 
axe  more  sensitive  to  changes  in  the  effective  electrical  conductivity  with 
temperature.  However,  they  are  essential  in  understanding  the  final  width 
of  the  spot. 

2.  Formulation. 

A  thin  cylindrical  ceramic  sample  is  positioned  in  a  single  mode  waveguide 
appUcator,  so  that  the  electric  field  along  its  axis  varies  as  the  fundamental 
mode  of  the  waveguide.  The  sample  is  held  in  place  at  its  ends  by  two 
thermally  insulated,  microwave  transparent  push-rods.  Although  the  electric 
field  is  altered  by  the  presence  of  the  ceramic,  for  the  present  analysis  we 
assume  that  this  effect  is  negligible  and  that  the  time-harmonic  electric  field 
is  given  by 

E  =  Eo  sin(7rZ/d)j  (1) 

where  d  is  the  height  of  both  the  cylinder  and  the  guide,  Eq  is  the  strength 
of  the  incident  mode,  and  j  is  a  unit  vector  perpendicular  to  the  axis  of 
the  cylinder.  This  assumption  effectively  decouples  the  equations  for  the 
electromagnetic  field  from  the  equation  for  the  energy  of  the  sample,  and 
allows  us  to  focus  solely  on  the  sample’s  thermal  field.  In  addition,  we  assume 
that  the  sample  is  thin  enough  to  ensure  that  variations  in  the  electric  field 
are  negligible  across  its  circulax  cross-section. 

In  light  of  these  assumptions,  the  temperature,  T,  satisfies  the  energy 
equation 

^ipCpT)  =  V-{KVT)  +  ^^E^sm^{nZ/d\  0  <  Z  <  d,  0<R<a 

at  2 

(2) 
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where  R  =  y/X"^  +  is  the  radial  distance  from  the  cylinder’s  axis,  p  is  the 
density  of  the  ceramic,  Cp  is  its  specific  heat,  K  is  its  thermal  conductivity, 
a  is  its  effective  electrical  conductivity,  and  a  is  the  radius  of  the  sample. 
Although  variations  of  the  thermal  parameters  K  and  C p  are  small  over  the 
temperature  range  required  for  sintering  or  joining  when  compared  to  the 
change  in  the  electrical  conductivity,  they  are  included  in  the  following  anal¬ 
ysis,  since  they  may  have  a  profound  effect  on  the  dynamics  of  the  heating 
process. 

We  also  require  that  the  temperature  sati.sfies  the  surface  heat  balance 

Ar^T-f-/i(T-T^)  +  se(T'' -T^)  =  0,  R  =  a,  Q<Z<d  (3a) 
or 

where  /i  is  a  constant  corresponding  to  heat  loss  from  the  surface  by  convec¬ 
tion,  5  is  a  constant  for  radiative  heat  loss,  e  is  the  emissivity  of  the  surface, 
and  Ta  is  the  ambient  temperature  of  the  surrounding  medium.  To  simplify 
the  analysis  that  follows,  we  assume  that  the  ambient  temperature  remains 
constant.  At  the  ends  of  the  sample,  we  prescribe  the  boundary  conditions 

d 

-^T  =  0,  Z  =  0, d;  and  0  <  i?  <  a  (3b) 

C/iO 

and  we  take  the  initial  temperature  of  the  sample  to  coincide  with  the  am¬ 
bient  temperature,  i.e., 

r(x,y,z,o)  =  r^.  (4) 

Equations  (2)-(4)  constitute  a  nonlinear  initial  boundary  value  problem  for 
the  temperature  T  within  the  sample.  The  nonlinear  character  of  this  sim¬ 
plified  problem  is  caused  by  the  dependence  of  the  electrical  conductivity  <t 
and  thermal  parameters  K  and  Cp  on  the  temperature  and  by  the  radiative 
losses  at  the  sample  boundary.  This  is  the  generalization  of  the  mathemat¬ 
ical  models  for  microwave  heating  as  studied  by  Tian  using  finite  difference 
simulations  [10]  and  by  Kriegsmann  using  asymptotic  methods  [5,6]. 
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3.  The  Simplified  Theory. 


There  axe  three  smaJl  parameters  that  arise  from  the  nondimensionaliza- 
tion  of  equations  (2-4).  The  first  is  the  Biot  number  Bi  =  ha/KA,  where 
Ka  =  K{Ta)  is  the  value  of  the  thermal  conductivity  at  the  ambient  temper¬ 
ature  Ta,  and  the  second  is  B2  =  seaT\l Ka-  The  former  is  a  measure  of  the 
relative  effects  of  convection  and  conduction  and  the  latter  is  a  measure  of 
the  relative  effects  of  radiation  and  conduction.  Typical  values  of  B\  and  B2 
for  ceramics  are  of  the  order  of  0.01,  see,  e.g.,  [7].  The  third  small  parameter 
is  the  fineness  or  aspect  ratio  of  the  cylinder  ajd. 

In  recent  studies  [3,4]  of  the  microwave  heating  of  ceramic  slabs  and  other 
compact  geometries  we  have  utilized  the  size  of  these  three  parameters  to 
obtain  an  asymptotic  approximation  to  the  temperature.  Similar  methods 
can  be  employed  for  the  nonlinear  initial  boundary  value  problem  (2-4)  but 
will  not  be  reported  in  detail  here.  The  result  of  this  analysis  is  that,  as 
Bi  — >  0  with  B2/B1  held  fixed  (so  that  the  effects  of  radiation  and  convection 
are  of  equal  importance)  and  with  (a/d)^ /B^  held  fixed,  the  temperature  of 
the  sample  T  is  given  by 

T{X,Y,Z,t)  =  U{Z,t)  +  0{B^)  (5) 

where  U{Z,t)  is  the  leading  order  approximation  to  the  temperature,  which 
is  independent  of  the  cross-sectional  or  transverse  coordinates  X  and  F,  and 
0{Bi)  represents  a  contribution  which  is  of  the  same  (small)  order  as  Bi. 

We  note  that  this  0(Bi )  error  term  does  depend  upon  X  and  F  and  thus 
represents  a  small  nonuniform  heating  across  the  cross-section  of  the  sample; 
the  term  will  not  be  calculated  here. 

The  leading  order  term  U{Z,t)  of  (5)  satisfies  the  dimensionless  initial 
boundary  value  problem 

—  (r(l-t-u))  =  — (fc— u)-|-p/(u).sin^(7r2)-2{u-|-Qr[(u-t-l)‘*-l]},  0  <  2  <  1, 

(6a) 

where  the  boundary  conditions  (3b)  and  initial  condition  (4)  become,  respec¬ 
tively, 

—  u  =  0,  z  =  0, 1  (6b) 
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lt{2,0)  =  0,  0  <  2  <  1.  (6c) 

In  arriving  at  equations  (6a-c)  we  have  introduced  the  dimensionless  variables 


2  =  Z/d, 


h 

a{pCp)A 


(7a) 


where  Z  is  nondimensionalized  with  respect  to  the  cylinder  length  d,  u  is  the 
relative  deviation  of  T  from  T/i,  and  t  is  nondimensionalized  with  respect 
to  the  ambient  convective  time  a{pCp)A/h.  We  have  also  introduced  the 
dimensionless  parameters 


cctaEq  B2  2 

and  the  dimensionless  functions 


fin) 


cxiiTAjl  +  «)) 

<^A 


k  =  K/Ka, 


pCp 

{pCp)A 


(7b) 


(7) 


where  a  a  =  <^{Ta)  is  the  effective  electrical  conductivity  at  the  ambient 
temperature,  a  is  the  ratio  of  the  convective  and  radiative  Biot  numbers, 
and  p  is  a  dimensionless  power. 

The  nonlinear  initial  boundary  value  problem  (6)  constitutes  the  mathe¬ 
matical  statement  of  our  small  Biot  number  theory  for  the  heating  of  the 
ceramic  rod. 


4.  Analysis. 

We  note  here  that  in  some  applications  where  ceramic  fibers  are  sintered 
in  a  single  mode  applicator  [11]  the  parameter  e  is  very  small.  For  other 
experiments  such  as  joining,  e  may  not  be  eis  small.  The  asymptotic  limit 
e  — »  0  is  quite  relevant  in  the  former  case  and  is  expected  to  give  qualitative 
results  in  the  later.  In  mathematical  terms,  the  theory  which  follows  is 
strictly  valid  for  the  ordering  Bi  <<  «  1.  That  is,  B\  «  a/d  «  y/B^, 

so  that  the  diffusion  term  in  (6c)  is  larger  than  the  neglected  terms  of  0{B\ ). 
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FIGURE  1 


Setting  c  =  0  in  (6a),  we  obtain  the  ordinary  differential  equation 

^(r(l  +  u))  =  p/(u)sin^(7rz)  -  2{u  +  a[(u  +  1)“*  -  1]}  (8) 

OT) 

the  solution  of  which  depends  upon  z  parametrically  and  satisfies  the  initial 
conditions  (6c). 

A  reasonable  model  for  the  effective  electrical  conductivity  leads  to  the 
function  /  being  given  by  the  Arrhenius-like  law  [12] 

/(u)  =  l-bc,e-‘^^/"  (9) 

where  Ci  ajid  C2  are  constants.  If  we  fix  z  and  define  P  =  psin^(7rz),  then  the 
solution  of  (8)  and  (6c)  increases  monotonically  from  its  initial  value  u  =  0 
to  a  terminal  value  u,  which  is  given  implicitly  by  the  solution  of 


P  =  G{v)  = 


2{t;-hft[(u-H)^  -  1]} 

fiv) 


(10) 
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A  graph  of  G(v}  is  shown  in  Figure  1,  from  which  we  deduce  that  there  can 
be  either  one  or  three  solutions  of  (10)  depending  upon  the  value  of  P.  If 
P  <  G,n  then  the  terminal  solution  v  lies  on  the  lower  branch,  whereas  if 
P  >  Gf4  then  it  lies  on  the  upper  branch.  If  G„,  <  P  <  G m  ,  then  there  are 
three  solutions;  one  on  the  upper  branch,  another  on  the  lower  branch,  and 
the  third  solution  on  the  middle  branch.  A  simple  analysis  of  (8)  shows  that 
solutions  on  the  upper  and  lower  branches  are  stable,  aaid  that  solutions  on 
the  middle  branch  are  unstable. 

We  observe  that,  because  of  the  spatial  variation  of  the  power  P  along  the 
axis  of  the  waveguide  and  sample,  at  different  points  along  its  axis  the  ce¬ 
ramic  sample  experiences  different  values  of  P,  and  so  there  is  the  possibility 
that  a  steady  temperature  distribution  may  be  on  the  upper  branch  in  one 
part  of  the  sample  while  it  is  on  the  lower  branch  in  the  remainder.  This  is 
indeed  the  case  if  we  tahe  the  dimensionless  power  p  >  Gm-  If  we  define  z\ 
by  ^ 

z\  =  —arcsin{y/GM/p)  and  Z2  =  I  —  zi,  (11) 

TT 

then  P  >  Gm  in  the  interval  z\  <  z  <  Z2,  and  a  steady  temperature  distri¬ 
bution  in  that  part  of  the  sample  must  lie  on  the  upper  branch.  We  have 
sketched  this  in  Figure  2. 

If  we  try  to  resolve  the  discontinuity  in  this  steady  state  approximation  by 
introducing  boundary  layers  at  zi  and  Z2 »  within  which  the  diffusion  term  of 
(6a)  is  important,  then  we  are  immediately  struck  with  the  disconcerting  fact 
that  such  a  solution  does  not  exist.  To  see  this,  we  introduce  the  stretched 

.  ,  ,  z  —  Zi  .  ,  ,  d  ,  ,  . 

variable  z  = -  into  (6a),  set  —  =  0,  and  obtain 

(  OTJ 

+  f{u)[GM  -  G{u)]  =  0,  Izj  <  cx)  (12a) 
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FIGURE  2 


where  G  is  defined  in  (10).  The  boundary  conditions  for  this  equation  are 
deduced  by  asymptotic  matching  of  its  solution  as  2  — >  ±oo  with  solutions  of 
(10).  We  find  from  a  str^lightfo^w^lrd  ajieilysis  that  u— +v„  +  o(l)as2— 
where  is  the  value  of  u  on  the  upper  branch  corresponding  to  Gm  (see 
Figure  1).  We  deduce  from  a  similar  einalysis  that  u  v\f  +  ^  as  2  — ^  —00. 
In  either  case  we  have 


— «  — ♦  0  as  2  — >  ±c». 
dz 


(12b) 


Upon  multiplying  (12a)  by  k  —  u,  integrating  the  result  from  ±00,  and  using 

dz 

(12b)  we  obtain 


fVu 

I  k{u)f{u)[GM—G{u)]du  =  0. 
JvM 


However,  from  the  definition  of  we  deduce  that  since  vm  <  u  <  Vu  the 
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integrand  is  positive  so  that  (13)  can  not  be  true,  and  a  steady  solution  of 
this  type  does  not  exist. 

We  resolve  this  apparent  contradiction  by  removing  the  constraint  that 
—  =  0  in  the  above  analysis.  That  is,  we  shall  look  for  a  traveling  wave 
solution  of  (6a)  which  has  the  form 


u  =  4>{z),  where  z  = 


and  T  =  erj. 


Inserting  this  ansatz  into  (6a)  we  obtain,  at  leading  order, 

^(*^*^)  +  fWlpsiii^  7r(zt  +C)-  G(<^)]  =  0  (15a) 

where  the  prime  on  C  denotes  a  derivative  with  respect  to  its  slow  time 
argument  t  =  erj,  aaid  the  equation  is  to  be  solved  on  the  interval  —  oo  < 
5  <  oo.  We  deduce  similar  boundary  conditions  to  the  above  and  find  again 
that,  as  in  (12b), 

d 

—  <f>  -^0  as  z  —y  ±oo.  (15b) 

dz 

Now  however,  asymptotic  matching  implies  that  the  solution  (p  has  the  limits 


(f>  —y  <t)±{C)  as 


(15c) 


where  </>+(C)  >  (p-{C)  are  those  roots  of  psin^(zi  +  C)  =  G{<f))  which  lie  on 
the  upper  and  lower  branch  of  the  S-shaped  response  curve  of  (10),  respec¬ 
tively.  Multiplying  (15a)  by  integrating  the  resulting  expression,  and 

dz 

applying  (15b)  and  (15c),  we  deduce  that 


’'(^1  +  c)  -  G(^))  d(^ 


ST^i‘W(:f/nr  +  t(i  +  ^)ldz 


where  the  dot  above  F  denotes  its  derivative  with  respect  to  (j). 
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Equation  (16)  is  a  first  order  nonlinear  ordinary  differential  equation  for 
the  position  of  the  slowly  moving  traveling  wave,  or  front,  which  ha-s  the 
initial  condition  C(0)  =  0  and  connects  the  solutions  on  the  upper  and  lower 
branches  of  (10).  We  note  that  the  dependence  of  the  right  hand  side  on  C  is 
implicit,  and,  in  particular,  the  integrand  in  the  denominator  depends  upon 
C  through  the  solution  of  the  boundary  value  problem  (15). 

We  can  now  deduce  the  dynamics  of  the  heating  process  when  p  >  G  ,  by 
using  (16).  Initially,  the  rod  heats  according  to  (8),  so  that  a  discontinuous 
temperature  profile  begins  to  form,  with  discontinuities  at  z  =  2]  and  z  =  Z2 
and  a  hot  spot,  where  z\  <  z  <  22,  as  shown  in  Figure  2.  At  this  point  in 
time,  considering  the  dynamics  of  the  left  half  of  the  sample  0  <  2  <  1/2 
alone  since  the  solution  is  symmetric  about  the  sample’s  midpoint  2  =  1/2, 
G  =  0  and  from  (16)  C  is  negative,  if  we  a.ssume  that  the  term  r+r(l  +  (^)  > 
0  as  will  be  the  case  for  physically  realistic  applications.  Thus,  C  decreases, 
and  the  front  begins  to  move  to  the  left,  from  z  —  z\.  This  elevates  a  larger 
portion  of  the  rod  to  the  higher  temperatures  of  the  upper  branch,  i.e.,  the 
hot  spot  begins  to  grow. 

This  description  is  given  under  the  proviso  that 

+ <«  >  0 

which,  from  the  definitions  of  and  F,  is  equivalent  to  the  statement  that 
the  internal  energy  density  of  the  ceramic,  pC pT,  is  an  increasing  function  of 
the  temperature.  This  is  true  in  ceramics  and  in  almost  eill  materials,  away 
from  phase  transitions. 

We  now  turn  to  the  further  development  and  stabilization  of  the  hot  spot. 
Noting  the  definitions  of  <A±(C),  it  follows  that  the  derivative  with  respect 
to  C  of  the  numerator  on  the  right  hand  side  of  (16),  i.e., 

r<t>+{0) 

N=  k{(t>)f{<f))\psiri^Tr{zi +C)-G{(f>)]d(l>,  (17) 

U-ic) 

—  iV  =  7rp/  fc((^)/(^)sin27r(zi  +  C)(iu,  (18) 
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which  is  strictly  positive,  since  0  <  zi  +  C  <  1/2,  and  k  and  /  are  positive. 
Also,  since  <^+(0)  =  </>-(0)  =  vm,  arid  =  psin^(7rzi ),  we  deduce  from 

(17)  that  N(0)  >  0.  Similarly,  if  we  define  C  so  that  G„,  =  psin^  7r(zj  +  C), 
then  ^+(G)  =  Vm,  <^-(G)  is  the  corresponding  point  on  the  lower  branch  of 
the  S-shaped  curve,  and  consequently  N{C)  <  0.  Hence,  there  is  a  unique 
C  =  C,  <  0  such  that  iV(C*)  =  0.  Since  the  denominator,  D(C),  of  (16)  is, 
for  the  reasons  given  above,  strictly  positive,  we  deduce  that  the  solution  of 
the  differential  equation  C  =  —N{C)/D{C)  with  initial  condition  C(0)  =  0 
is  monotone  decreasing  and  tends  to  C»  as  t  ^  oo,  with  its  final  approach 
being  exponential  in  r.  Thus,  the  hot  spot  grows  in  size  and  finally  stabilizes 
to  occupy  the  region  21  +  C,  <  2  <  1  —  —  C*. 

We  can  now  consider  the  influence  of  a  temperature-dependent  thermal 
conductivity  on  the  equilibrated  value  C*.  First  we  recall  that  G,  is  such 
that 


N{Ci,)=  I  k(<t>)f(<f>)[psin^  n{2i  +  C^)  -  G{(t>)]d(l>  =  0,  (19) 

and  consider  the  case  when  the  thermal  conductivity  is  constant,  so  that 
k(4i)  —  1.  Then,  since  and  are  such  that  the  local  power  at 

the  boundaries  of  the  hot  spot  P  =  psin^(2i  -fC*)  =  G(<^±(C*)),  the  graphs 
of  /(it)G(u)  and  Pf{u)  intersect  transversally  at  u  =  <^_(C*),  u  =  <^+(C*), 
and,  from  (19),  also  at  some  value  between,  i.e.,  on  the  interval  <^_(C,)  < 
u  <  <^+(C*).  Equation  (19)  implies  that  the  unique  value  C*  is  such  that 
the  areas  of  the  two  lobes  between  the  graphs  of  /(u)G(u)  and  Pf{u)  are 
equal.  When  k{4>)  is  not  constant  but,  for  exeimple,  is  a  monotone  increasing 
function  of  (f)  with  k{0)  =  1,  as  is  the  case  for  typical  ceramics,  the  influence 
of  the  temperature-dependence  with  C  fixed  is  such  as  to  increase  the  area 
of  the  right  hand  lobe  (at  larger  <^)  more  than  that  of  the  left  hand  lobe. 

This  increases  iV(C),  so  that,  since  >  0,  the  temperature-dependence 

of  k  is  such  as  to  decrease  the  equilibrating  value  C*  to  more  negative  values, 
and  hence  increase  the  final  width  of  the  hot  spot  in  the  steady  state. 
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5.  Conclusion. 


The  implications  of  the  analysis  for  the  sintering  of  ceramic  fibers  and 
joining  of  ceramic  cylinders  is  now  evident.  In  the  first  case,  the  hot  spot 
forms,  propagates,  and  then  stabilizes.  If  the  temperature  in  the  relatively 
warm  region  of  the  hot  spot  is  sufficient  for  sintering,  then  the  fiber  can 
be  slowly  pulled  through  the  guide,  thus  insuring  that  the  entire  sample  is 
processed.  This  is  used  as  a  means  of  sintering  ceramics  in  practise  [11],  and 
the  rate  at  which  the  fiber  is  to  be  drawn  is  found  experimentally.  In  the 
second  case,  the  hot  spot  is  to  encompass  the  butt  joint  at  which  the  two 
ceramic  cylinders  are  to  be  joined,  and,  if  the  temperature  in  this  region  is 
sufficient  for  the  materials  to  fuse,  then  a  strong  joint  can  be  obtained  [2]. 

We  close  by  briefly  describing  another  type  of  solution  that  is  possible  if 
the  applied  electric  field  has  a  minimum  at  the  center  of  the  fiber.  This  may 
occur  by  exciting  the  applicator  in  one  of  its  higher  spatial  modes.  If  the 
maximum  of  the  electric  field  is  such  that  P  >  Gm  and  the  minimum  is  such 
that  P  <  Gm,  then  hot  spots  will  form  at  both  eirds  of  the  fiber.  These 
spots  will  grow  in  size  and  stabilize  according  to  the  mechanism  described 
at  the  end  of  the  last  section. 
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